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SUMMARY 

Through the use of equivalent variable analogies, the authors demonstrate how an electrical subsystem can 
be modeled by an equivalent structural subsystem. This allows the electrical subsystem to be probabilistically 
analyzed by using available structural reliability computer codes such as NESSUS. With the ability to analyze 
the electrical subsystem probabilistically, we can evaluate the reliability of systems that include both structural 
and electrical subsystems. Common examples of such systems are a structural subsystem integrated with a health- 
monitoring subsystem, and smart structures. Since these systems have electrical subsystems that directly affect 
the operation of the overall system, probabilistically analyzing them could lead to improved reliability and reduced 
costs. The direct effect of the electrical subsystem on the structural subsystem is of secondary order and is not 
considered in the scope of this work. 


INTRODUCTION 

Complex systems are often composed of subsystems of more than one engineering discipline (e.g., electrical, 
mechanical, chemical, nuclear, fluid, and stmctural). Although each of the subsystems can affect the total system 
and often affect other subsystems, the components that comprise the subsystems may vary. For example, a struc- 
tural subsystem has material property, load, and dimension components, whereas a simple electrical subsystem 
has resistance, voltage, and current components. Each component in each of these subsystems varies in magnitude 
and/or direction. Because of the variability of the components and the various levels of sensitivity of their coupling, 
the complete system varies in its operations. In order to analyze the variability of such a complex system, the 
subsystems must be evaluated probabilistically (stochastically). And probabilistic evaluation of multidisciplinary 
systems requires computer codes that can analyze structural, thermal, chemical, fluid, and electrical subsystems. 
However, it is through such an evaluation that critical subsystems and components can be identified; this, in 
turn, can lead to reduced cost and improved reliability of the complex systems. 

One probabilistic code used to evaluate structural subsystems is NESSUS (Numerical Evaluation of Stochas- 
tic Structures Under Stress); it has been developed for the NASA Lewis Research Center by Southwest Research 
Institute (refs. I and 2). This code can be used to evaluate structural subsystems with variable material properties, 
dimensions, and loadings. Possible modes of failure can be combined through fault tree analysis to determine 
the reliability (probability of failure) of the subsystem. Risk analysis that includes the subsystems’ costs can also 
be performed with the code. However, this probabilistic structural code is limited to the analysis of structural 
and thermal-stmctural subsystems. 

The work reported herein draws analogies between the variables of electrical and structural subsystems and 
then uses them to analyze several simple electrical circuits. These equivalent variables allowed an electrical 



subsystem to be probabilistically analyzed with the NESSUS computer code. This ability to analyze the sub- 
system probabilistically allowed us to evaluate the reliability of systems having both structural and electrical 
subsystems. In the future, identification of the critical components, whether structural or electrical, will make 
possible more reliable systems at reduced costs. 

Common examples of such systems are a structural subsystem integrated with a health-monitoring subsystem 
(ref. 3), and smart structures (ref. 4). These systems have electrical subsystems that directly affect the operation 
of the overall system. For example, in a health-monitoring subsystem, electrical sensors inform system operators 
to make proper adjustments. If the monitoring system fails, the structural subsystem can be damaged or it can fail. 
Similarly, in smart structures that have an electrical subsystem in the form of piezoelectrical sensors and/or actuators 
(ref. 5), failure of the piezoelectric system can also result in failure of or damage to the structural subsystem. 

The primary objective of this work was to determine the feasibility of using a structural code to probabilis- 
tically analyze electrical subsystems. Investigating the direct effect of the electrical subsystem on the structural 
subsystem, however, was beyond the scope of this work. 

Analogies between other disciplines and the structures discipline can be demonstrated so that such subsystems 
can be probabilistically analyzed with a single code. 


ELECTRICAL AND STRUCTURAL ANALOGIES 

There are analogies between the variables in the equations used by various engineering disciplines; thus, 
similar equations can be used amongst the disciplines. Analog computers were the prime example of the use of 
such analogies (ref. 6). Here, in an inverse example, dc and ac circuits are analyzed by using structural analogies 
and computer codes developed to analyze structural systems. 

A dc circuit can be analyzed by the following equation: 


R 


where I is the current, V the voltage, and R the resistance. The equivalent equation for structured subsystems is 


F = ku 


( 2 ) 


where F is the force, k the stiffness of the structure, and u the displacement. Equations (1) and (2) show that 
force F is analogous to current /; displacement u is analogous to voltage V; and flexibility l/k is analogous to 
resistance R. Though other combinations of analogies can be derived for the variables in these equations, those 
just presented result in consistent analogies, which will be used here. In addition, these analogies have a physi- 
cal correlation in that force flows in a structure just as current flows in a circuit; the portion of force transferred 
to members depends on their flexibility just as the portion of current transferred to various paths depends on 
their resistance. And displacement is a function of force and flexibility just as a voltage change is a function of 
current and resistance of the circuit. Table I summarizes the variables that are equivalent in the simple electrical 
and the static structural subsystems. The makeup of these electrical circuits is limited to voltage and current 
sources and resistors. 


2 



For electrical subsystems that include inductors and capacitors, a dynamic structural subsystem (ref. 7) can 
be used to model the electrical subsystem. The relationship for an inductor is 


I = 



( 3 ) 


where L is the inductance, and t is time. The similar equation for the structural model is 

= ku = k j li dt 

where it is the velocity and F, is the stiffness force. Comparison of equations (3) and (4) shows flexibility is 
equivalent to inductance L, and velocity u is analogous to voltage V. 

The equation for a capacitor is 



(5) 


where C is the capacitance. The analogous equation for a structure is 

Fj = -mit 


where F, is the inertial force, and u is the acceleration. In equations (5) and (6), mass m is analogous to 
capacitance C. 

The structural equation that is equivalent to the resistor equation in such a circuit is 

F, = ai (7) 


where F^ is the damping force, and c is the damping. From equations (1) and (7) we can see that c is analogous 
to the inverse of the resistance. 

The linear equation to be solved in the electrical circuit is 


dt R L J 


( 8 ) 


which is analogous to the structural dynamics equation 

mu+ cii + ku = F{t) 

where F is the forcing function. The equivalent variables for the ac circuit are summarized in table II. 
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In both of the foregoing analogies, the force is taken to be analogous to the current. However, force is a 
vector and current is a scalar. Therefore, uniaxial structural models were used to eliminate the difficulties asso- 
ciated with the direction of forces when an electrical circuit was modeled. It is possible to use planar structures 
or nonplanar structures to model more complex electrical circuits. When nonplanar structures are used, care 
must be exercised to properly represent the analogies by the structural member’s spatial position. The use of 
nonplanar electrical circuits, however, introduces comparable complexities in the analogy process. 


Direct Current Circuit Analyses 

The simple electrical circuit shown in figure 1(a) was probabilistically analyzed with the NESSUS computer 
code. The circuit was modeled as an equivalent structural subsystem by using the analogies previously presented. 
Figure 1(b) shows the structural model with node and element (encircled) numbers. Elements 1, 4, and 5 had 
nondimensional flexibilities of 0.001, 0.004, and 0.005, respectively. These flexibilities correspond to the resis- 
tors of the circuit. Elements 2 and 3 were made 10(X) times stiffer than element 1 to model a rigid member and 
allow the proper flow of forces (current) to members 4 and 5. All elements were beam elements. Therefore, 
only axial degrees of freedom were allowed in order to prevent bending. Node sets (2,3), (4,6), and (5,7) are 
master and slave node sets for the different member sections connected at these locations. The model was 
subjected to a normalized displacement of 12 at nodes 8 and 9 to simulate the voltage source. 

Several analyses were performed with the model. The stiffness of the elements (equivalent resistance) was 
varied by perturbing the cross-sectional height and the modulus of elasticity of the elements. Although either of 
these two variables could have been perturbed to model the variability in the stiffness (resistance), both were 
perturbed to determine their combined effect on the results. The assumed distributions for the perturbed variables 
are shown in table III along with assumed means and scatter expressed in percent of the mean. The scatter of 
the variables signifies expected fluctuations of the variable. 

The cumulative probabilities of occurrence for the currents and for the voltage drop V, are shown in fig- 
ures 2 and 3. The 95th percentiles show values of 3763, 2115, and 1693 A for currents I, I 2 , and I 3 , respectively, 
and 3.79 V for voltage V,. These are values that have only a 5-percent chance of being exceeded. Deterministic 
analysis resulted in values of 3724, 1655, and 2069 A and 3.72 V for currents I,, I 2 , 13, and voltage V,, respec- 
tively. The deterministic veilues are approximately at the 50th percentile level for the cumulative distribution 
functions. From this comparison, we observe that slight differences will occur if the design of the circuit is 
based on a deterministic analysis rather than on the 95th percentile of a probabilistic analysis. However, a deter- 
ministic analysis does not provide the designer with information on the probabilities that the responses will be 
exceeded, and hence, on the reliability of the design. 

Sensitivities are a measure of the effect that each variable has on the probability of the response. Averages 
of the sensitivities for the perturbed variables are depicted in figures 4 and 5. The averages were used because 
the sensitivities varied little at various levels of probability. 

The cross-sectional height and modulus of elasticity of the structural elements correspond to the flexibility, 
which is analogous to the resistance within the circuit. This analogy allows the notation (see figs. 4 and 5) for 
referring to the resistor perturbation and to the corresponding structural perturbations. For example, the sensitiv- 
ity for resistance 1 (Rl) caused by perturbing the analogous element’s cross-sectional height is resistance 1-h 
(Rl-h). Similarly, the Rl sensitivity induced by perturbing the modulus of elasticity of the element is represented 
by resistance 1-m (Rl-m). Therefore, both resistance 1-h and resistance 1-m should be considered in drawing 
conclusions about the sensitivity in the analysis caused by resistance 1 . 
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Figure 4 shows that in the analysis for I, the sensitivities of the variables related to R1 have a considerable 
effect. The height variable analogs of resistors 2 and 3 also have a significant effect on I,, but the modulus of 
elasticity analogs do not. It can also be concluded from this figure that R3 has a greater effect on the I, analysis 
that R2 does. In the analysis for Ij, the sensitivities of the variables show that R2 has the greatest effect on Ij, 
followed by Rl. As would be expected, R3 has little effect on I 2 . Sensitivities of the variables in the analysis for 
I 3 show just the opposite results, with R3 having the greatest effect, Rl the next greatest, and R2 the least. In the 
analysis of V„ sensitivities show results similar to that of I„ but with even more of an effect coming from Rl. 

The other voltages can be determined by examining the relative displacements (between nodes 7 and 9 for 
element 4, and nodes 6 and 8 for element 5) as opposed to absolute displacements. Although determining these 
voltages was not within the scope of this study, there is no fundamental limitation of the codes that would pre- 
vent such a determination. 

A separate analysis was done to determine the probability of the voltage V, exceeding an allowable voltage 
(equivalent displacement) that was assumed to have a scatter range described by the distribution and parameters 
shown in table III. This resulted in a probability of failure of 0.18482. If the allowable voltage were considered 
deterministic, with a value of 3.9028 V, the probability of failure would drop to 0.0000019. This is shown in 
figure 3 by the high probability that the resultant voltage would be less than the allowable voltage. Thus, includ- 
ing the scatter of the allowable significantly affects the analysis. 

A circuit with a current source was also probabilistically analyzed. It and the equivalent structural model are 
shown in figure 6 . Elements 1, 2, 5, and 6 were 500 times stiffer than the next stiffest element, 8 . Again, this 
was done to make these elements as rigid as possible, thereby allowing proper transfer of the forces (currents). 
Elements 3 and 4 correspond to 10-f2 resistors, whereas elements 7 and 8 correspond to 20- and 5-f2 resistors, 
respectively. A concentrated nondimensional load of 36 was placed at node 1 to simulate the current source. All 
nodes were allowed only horizontal translation except nodes 7, 12, and 14, which were fixed. These boundary 
conditions removed bending effects to accurately model the circuit. 

Several analyses were performed on this model. Table IV shows the distributional information about the per- 
turbed variables used in the analyses. The cumulative probabilities of occurrence for currents I,, I2, I3, and I4 are 
shown in figure 7, and the average sensitivity factors are expressed in figure 8 . The 95th percentile values for 
currents I„ I 2 , 13, and I4 were 24.47, 17.47, 13.98, and 3.51 A, respectively. Deterministic analyses resulted in 
corresponding values of 21, 15, 12, and 3. The difference between the probabilistic and deterministic analyses is 
approximately 16.5 percent, but the deterministic analyses provide no information on the chance that the responses 
will be exceeded. 

The sensitivities of the currents in the analyses showed that the current source had a very significant effect 
on all analyses of the currents. The resistors in which the current was found had the next largest effect. Other 
variables had an insignificant effect. 

The cumulative probabilities of occurrence of voltages V„ V3, and V4 are shown in figure 9, and average 
sensitivity factors are shown in figure 10. Voltage V 2 was not determined for reasons previously discussed. The 
curves for Vj and V 4 overlap and are nearly identical, as would be expected. The 95th percentile voltages for 
V|, V 3 , and V 4 were 245.25, 69.83, and 70.08, respectively. Deterministic analyses resulted in values of 210, 60, 
and 60 for the same voltages. This is approximately a 16.8-percent difference between the probabilistic and 
deterministic analyses. 

In this analysis, the sensitivities of the voltages showed that the current source had an even more significant 
effect than it did in the analyses of the currents. The sensitivities of all the resistors had an almost negligible 
effect. 
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Alternating Current Circuit Analysis 


A simple circuit (fig. 11(a)) containing a resistor, inductor, and capacitor was probabilistically analyzed. This 
was done by modeling the circuit as the structural system shown in figure 11(b) and using the NESSUS com- 
puter code. Because of the presence of a capacitor and inductor, a harmonic analysis was performed by using 
the analogies listed in table II. A harmonic forcing function with an amplitude (A) of 0.667, a phase angle ((|)) of 
idl rad, and a frequency (to) of 1.0 rad was placed at node 1 to simulate the current source. The governing 
equation for this problem is 

mu + cu + ku = 


By choosing these forcing parameters, we modeled the imaginary current source shown in figure 11(a). Two 
beam elements were used to improve the solution even though only one element is required to model all three 
of the electrical components. To prevent any bending effects, only horizontal translation of nodes 1 and 2 was 
allowed. A lumped mass of 0.0001 was placed at node 1 to model the capacitor. This required that the mass 
density of the elements be very small so that the results would not be adversely affected. The total axial stiff- 
ness of the structure was set at 0.0333 to model the inductor. A damping factor (Q of 0.43651 was used; it 
relates to the resistor through equation (7) as follows: 



( 11 ) 


where c is the damping (1/resistance), k the stiffness, and m the mass. The variables that were perturbed in the 
model are shown in table V with their corresponding distribution information. Inductance- 1 is analogous to 
perturbing the length of the structural members, which is its flexibility and is analogous to inductance. 

The structure was analyzed twice. The first analysis determined the cumulative probability of occurrence for 
the imaginary component of the voltage (velocity). The second analysis determined the probability that the real 
component of the voltage would occur. The cumulative probability of occurrence results and the results of the 
sensitivity analyses are shown in figures 12 to 14. Probabilistic analyses resulted in 95th percentile voltages of 
1.02 and 20.90 for the imaginary and real components of the voltage, whereas deterministic analyses resulted in 
values of 0.96034 and 20.024, respectively. There is approximately a 5-percent difference between the values 
obtained with the two types of analysis. 

The sensitivities show that the variables related to inductance had the greatest effect on both analyses, and 
the current source had the next greatest effect. The resistance affected the analysis for the imaginary component 
of the velocity slightly, but had no effect on the analysis for the real component. 


POTENTIAL ENHANCEMENTS 

Enhancements to the structural computer code and/or modeling procedure would allow detailed analyses of 
complex electrical circuits. In this investigation, an effort was made to identify possible enhancements without 
incorporating them in the code. Some are described here. 

Modifying the code to allow boundary conditions to be perturbed would permit analysis of the effect of a 
stochastic voltage source that would be modeled with a prescribed boundary displacement. Using springs as supports 
that can be perturbed might be a modeling remedy. However, the correct spring stiffness is not obvious. 
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Models for overdamped electrical circuits (i.e., damping factor >1) are not accounted for by comparable 
structural analogs, and appropriate enhancements to the structural code will have to be made in order to model 
overdamped electrical circuits. 

An enhancement that permits the analog for conductance, with an additional member that has no stiffness 
but has a density, is needed. 

The direct input of the damping matrix may be a solution to analogous structures that do not have a constant 
damping ratio. This would apply to structural analogs for modeling parallel circuits (fig. 11). 

Slight modifications to the code to allow the inertia, damping, and stiffness forces to be valid response out- 
puts would also permit probabilistic analyses of the currents within the conductor, resistor, and inductor. 

These enhancements would make it possible to probabilistically analyze structural systems that contain elec- 
trical systems such as health-monitoring systems and smart structures. Additional enhancements could be made 
to include the coupling effect between the electrical and structural systems. This would improve probabilistic 
system analyses and would lead to improved designs, from a safety as well as an economic viewpoint. 


CONCLUDING REMARKS 

We have shown that an electrical circuit can be probabilistically analyzed by using equivalent variables in 
the structural analysis code NESSUS. Therefore, system reliability analyses can be performed on systems con- 
taining both structural and electrical subsystems. There are other analogies that exist between engineering disci- 
plines that should allow other types of subsystems to be probabilistically analyzed with the same code. The end 
result would be system reliability analyses of a complete system, including such subsystems as electrical, thermal, 
hydraulic, structural, and the like. Analyses of coupling between subsystems could eventually be included, thereby 
leading to realistic modeling of the assembled system. Such improved system analysis would lead to increased 
reliability and reduced cost of the system. 
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TABLE L— ANALCX30US VARI- 
ABLES IN DIRECT CURRENT 
CIRCUITS AND STATIC 
STRUCTURES 


Electrical 

variables 

Equivalent 
structural variables 

Current, I 
Voltage, V 
Resistance, R 

Force, F 
Displacement, u 
Flexibility, Uk 


TABLE n.— ANALOGOUS VARIABLES 
IN ALTERNATING CURRENT 
CIRCUITS AND DYNAMIC 
STRUCTURES 


Electrical 

variables 

Equivalent structural 
variables 

Current, / 
Voltage, V 
Resistance, R 
Capacitance, C 
Inductance, L 

Force, F 
Velocity, U 
1/Damping, I/c 
Mass, m 
Flexibility, Uk 


TABLE m.— RANDOM VARIABLE STATISTICS 


FOR DIRECT CURRENT CIRCUIT 1 


Variable 

Distribution 

Mean 

Scatter, 

percent of mean 

Resistance 





1-h 

Normal 

0.5 

10.0 

2-h 

Normal 

.5 



3-h 

Normal 

.5 



l-m 

Lognormal 

2000.00 



2-m 

Lognormal 

400.00 



3-m 

Lognormal 

500.00 

' 


Voltage 1* 





(allowable) 

Normal 

3.9028 

10.0 1 


“Used in one analysis only. 
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TABLE IV.— RANDOM VARIABLE STATISTICS 


FOR DIRECT CURRENT CIRCUIT 2 


Variable 

Distribution 

Mean 

Scatter, 

percent of mean 

Current source 

Normal 

36.0 

10.0 

Resistance 







1-h 

Normal 

0.5 

10.0 

2-h 







3-h 







4-h 







l-m 

Lognormal 

.2 



2-m 



.2 



3-m 



.4 



4-m 




1 




TABLE V.— RANDOM VARIABLE STATISTICS 
FOR ALTERNATING CURRENT CIRCUIT I 


Variable 

Distribution 

Mean 

Scatter, 

percent of mean 

Current source 

Lognormal 

0.667 

10,0 

Resistance 

Lognormal 

.4365 

5.0 

Inductance-h 

Normal 

.5 

10.0 

Inductance-I 

Lognormal 

2.0 

1.0 

Inductance-m 

Lognormal 

.1333 

10.0 
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Figure 1. Direct current circuit 1 and analogous structural 
model, (a) Electrical dc circuit 1. (b) Structural analog 
showing node and element (encircled) numbers. 
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Figure 2. Probabilities of occurrence for currents of dc 
circuit ^ . 
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Figure 3. Probability of occurrence for of dc circuit 1 . 
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Figure 4. Sensitivities (relative units) for currents of dc circuit 1. 
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0.001 n 



Resistor length Resistor diameter 

Figure 5. Sensitivities (relative units) for of do circuit 1. 



(b) 

Figure 6. Direct current circuit 2 and analogous structural 
model, (a) Electrical dc circuit 2. (b) Structural analog 
showing node and element (encircled) numbers. 
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Resistor length Resistor diameter 


Figure 8, Sensitivities (relative units) for currents of dc circuit 2. 
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Figure 9. Probabilities of occurrence for voltages of dc 
circuit 2. 
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Figure 10. Sensitivities (relative units) for voltages of dc 
circuit 2. 
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(a) 


j— Lumped mass = 0.0001 



1 2 3 
(b) 


Figure 11 . Alternating current circuit and analogous structural 
model, (a) Electrical ac circuit, (b) Structural anolog. 
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Figure 14. Sensitivities (relative units) for voltages of ac circuit. 
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